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ABSTRACT 



o 
u 

■ We present polarization measurements at 8.4, 22, and 43 GHz made with the VLA 

of a complete sample of extragalactic sources stronger than 1 Jy in the 5-year WMAP 



Q ■ catalogue and with declinations north of —34°. The observations were motivated by the 

^ I need to know the polarization properties of radio sources at frequencies of tens of GHz 

^ ■ in order to subtract polarized foregrounds for future sensitive Cosmic Microwave Back- 

ground (CMB) experiments. The total intensity and polarization measurements are 
generally consistent with comparable VLA calibration measurements for less- variable 
sources, and within a similar range to WMAP fluxes for unresolved sources. A further 
paper will present correlations between measured parameters and derive implications 
' for CMB measurements. 

■ Key words: radio:polarization 

ON 

O 1 INTRODUCTION 



Radio source samples, selected at high frequency and with flux densities of the order of 1 Jy, mostly contain compact 
flat spectrum objects associated with the nuclei of active galaxies which have relativistic outflows pointing close to 



X 

d \ the direction of the observer (e.g. Blandford & Rees 1978; Scheuer & Readhead 1979; Orr & Browne 1982). This 
emission consists of a number of synchrotron components, each with a steep optically thin synchrotron spectrum at 
high frequencies and a self-absorbed spectrum at lower frequencies; the combination of many such components with 
different cutoff frequencies produces the observed "flat" radio spectrum. These components are frequently linearly 
polarized, and in principle useful information on jet physics can be obtained from the distribution of polarization in 
frequency and space. Although high frequency-selected samples are dominated by flat spectrum objects the samples 
do contain a minority of intrinsically more powerful sources in which the extended, steep-spectrum synchrotron 
emission dominates rather than the emission from a beamed core. Such sources also display linearly polarized 
emission, often at a higher level than the flat spectrum sources. 

In recent years there has been particular interest in high-frequency properties of radio sources because discrete 
radio sources act as a confusing foreground for measurements of the Cosmic Microwave Background (e.g. Tegmark & 
Efstathiou 1996). Such sources contribute most to the power spectrum on small angular scales (high £), and efforts 
have been made to subtract their effects in several different ways. For total intensity, high frequency observations of 
sources selected from lower frequency surveys have been used. For example, the Very Small Array CMB interferom- 
eter, working at 31 GHz, used 15-GHz source-finding observations with the Ryle Telescope followed by monitoring 
of discrete sources with dedicated long baselines (Dickinson et al. 2004) and/or single dish observations made with 
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the Torun 32 m (Gawronski et al. 2009), while for the Cosmic Background Imager experiment high frequency ob- 
servations were made with the OVRO 42 m and GBT of sources selected at 1.4 GHz (Mason et al. 2003; Mason 
et al. 2009). Alternatively, population models based on extrapolations from lower frequency catalogues can be used 
statistically to ameliorate the effects of the extragalactic source foreground. 

Although large-scale blind surveys at high frequencies are difficult to do, because of the small field of view at 
high frequencies, there arc a number of both blind and pointed surveys in existence. The southern hemisphere has 
been completely surveyed by the AT20G survey (Ricci ct al. 2004), which has been conducted with the Australia 
Telescope Compact Array (ATCA) at 18 GHz down to a limiting flux density of 50 mJy. A detailed study of brighter 
sources in this survey has been presented by Massardi et al. (2008). The whole sky has been covered to a brighter 
limiting flux density of about IJy by WMAP at a number of frequencies >20 GHz (Wright et al. 2009), and smaller 
areas of the northern sky have been covered to greater depth by surveys such as the 15-GHz 9C survey (Waldram 
et al. 2003). These surveys generally find that sources that are bright at high frequencies have complex spectra. 
About 30% arc significantly inverted (Ricci ct al. 2004), and the mean spectral index around 20-40 GHz tends to be 
around zero, generally steepening at higher frequencies (Ricci et al. 2006; Massardi et al. 2008; Vol'Vach et al. 2008). 
There is some evidence that the spectra of most sources turn over by 95 GHz (Sadler et al. 2008) but in general 
extrapolation to both higher frequencies and to lower flux density levels is needed in order to realistically simulate 
the CMB discrete foregrounds (e.g. Waldram et al. 2007). Importantly, there is convincing evidence that the mix of 
flat and steep spectrum sources changes with flux density. This is a prediction of beaming models (e.g. Orr & Browne 
1982; Wall & Jackson 1997). There is also direct observational evidence provided by Gawronski ct al. (2009) that the 
proportion of steep-spectrum sources with extended radio emission increases from <20% in a WMAP 22-GHz sample 
to >50% in Ryle-selected 15-GHz sample which has a flux density limit about two orders of magnitude weaker than 
that of the WMAP sample. Thus care must be taken when extrapolating to both higher frequencies and lower flux 
densities. The sensitivity of instruments such as Planck require this to be done as accurately as possible, and hence 
it is important to base the models on as firm an observational footing as possible. 

The situation with polarized foregrounds is less developed but it is becoming urgent to solve because inflation 
models make different predictions for the strength of B-modes in the CMB polarization distribution. A heroic effort is 
going into experiments to detect B-modes. These include Planck and ground-based projects like QUIET (Samtleben 
2008) and the balloon-borne experiment SPIDER (Montroy et al. 2006). All CMB observations are technically 
challenging because they rely on extreme stability of the equipment and on meticulous subtraction of foregrounds 
to get at the true distribution of the CMB fluctuations. The detection of the B-mode signal will be particularly 
difficult (though just how difficult depends on its strength which varies over a wide range depending on the inflation 
model). For scalar-to-tcnsor ratio r=0.01, the peak signal (at around 1° scales) is ~30nK. B-mode experiments 
operate at high frequencies (e.g. QUIET: 44 and 90GHz; Spider: 90, 145 and 220 GHz; PLANCK 30,70, 100, 143, 
217 and 353 GHz). Hence observations at the highest frequencies possible are desirable. Estimates of synchrotron 
source total intensity contamination have been made by extrapolating from observations made at 1.4 and 5 GHz (e.g. 
Tucci et al. 2004). The current best-guess models are presented by Toffolatti et al. (1998, 2005). However, additional 
assimiptions need to be made if the polarized foreground is to be characterized; one needs to know about the high 
frequency polarization properties of sources. Ideally one would like direct observations of polarized source counts but 
this is impractical. An alternative is to use high frequency total intensity source counts and a knowledge of how the 
percentage polarizations of sources change with both frequency and flux density to predict polarized source counts. 

The work reported in this paper is the first phase in a programme aimed at improving knowledge about the 
polarization properties of radio sources at high radio frequencies. We begin with bright sources from the WMAP 
sample of Wright et al. (2009) which have 6'22GHz >1 Jy- In this work, all but four of the sample of 203 WMAP 
sources with flux density of >l Jy at 22 GHz and with declination 5 > —34 have been imaged with the VLA in 
total intensity and polarization. The second phase of our programme will be to measure the polarization properties 
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of a high frequency-selected sample of much weaker sources designed to look for any dependence of polarization 
properties on flux density. Given that Gawronski et al. (2009) find that the population mix changes with flux density 
there is every reason to expect that average polarization properties also change. 

2 OBSERVATIONS AND DATA REDUCTION 

Observations were conducted using the VLA in D-conflguration between 01.40 and 19.30 UT on 2008 August 1. 
During this period 16 antennas were equipped with EVLA receivers while the remaining 11 had original VLA 
receivers. Each source was normally observed at three frequencies: two 50-MHz bands centred on 8.4351 and 8.4851 
GHz (X-band), two 5G-MHz bands centred on 22.4351 and 22.4851 GHz (K-band), and two 50-MHz bands centred 
on 43.3149 and 43.3649 MHz (Q-band). Pointing calibration was carried out approximately every 15-20 minutes. 
The available time allowed observations to be made of all the sources at 22 GHz and 43 GHz, and approximately 
two-thirds of the sources at 8 GHz. Four sources (WMAP 0823-^224, WMAP 0824+392, WMAP 1014+231 and 
WMAP 1840+797) were missed during the observations, leaving a total of 199 sources out of 203 with attempted 
measurements. The source 3C48 was observed in order to provide an additional flux and polarization calibrator; the 
primary flux and polarization calibrator, 3C286 (WMAP 1331+305) was already part of the sample. Additionally, 
three much weaker sources from the 9C survey were also observed as a pilot for a future observing programme aimed 
at much weaker sources. 

The nominal time on each source was 37 seconds at 8 GHz, 30 seconds at 22 GHz and 80 seconds at 43 GHz, 
giving a theoretical rms of 0.2, 0.5 and 0.7 mJy/beam, respectively, for natural weighting. However, the rms noise in 
most images is typically greater than this, being between 1 mJy/beam and 2 mJy/beam in the majority of images. 
This is partly due to the use of uniform weighting during the imaging process, which increases the noise level typically 
by factors between 1.2 and 1.5. However, in particular, 8-GHz observations were affected by a combination of missing 
antennas and slightly longer than predicted slew and setup times on many sources. A sampling time of 3 seconds 
per integration was used, the minimum allowed by the observing system. 

All data analysis and processing was carried out in the aips package, distributed by the US National Radio 
Astronomy Observatory. Data were examined and flagged manually. The 22-GHz data were relatively clean, but 
extensive flagging of bad telescopes was performed source by source on both 8-GHz and 43-GHz data, resulting in 
the removal of up to half of the data on any given source. The major problems were loss of correlation on some 
telescopes during setup at the beginning of a scan, together with intermittent bad data on many individual telescopes. 

Overall amplitude calibration was performed using observations of the source 3C286 at each frequency and 
normalising to the flux scale of Baars et al. (1977). Amplitude and phase solutions were made using the sources 
stronger than a limiting flux density; the process was iterated, adjusting this limiting flux, until the amplitude and 
phase solutions appeared to vary in a coherent manner. Amplitude solutions were smoothed using an averaging 
time of 1 hour for 43 GHz and 0.1 hours for the other frequencies. In the case of the 8-GHz observations, a further 
calibration was made, again using a standard model for 3C286, for baseline-dependent offsets. This was necessary 
because of the mixture of receivers in the array, approximately equally distributed between old (VLA) and new 
(EVLA) receivers. 

Polarization calibration was then carried out. Telescope polarization offsets were calibrated by use of the source 
3C84, which was assumed to have zero polarization. This assumption is correct for the level of accuracy required 
here: Taylor et al. (2006) report levels of <0.1% polarization in the core of the source at frequencies less than 
22 GHz, with approximately 0.2% at 22GHz, and a somewhat higher level (0.8-7.5%) in the weak mas-scale jet 
component. The overall level of polarization in 3C84 is likely to be about 0.1%. 3C84 was not observed at 8 GHz. 
In this case the source 0125—001, a bright point source for which no signiflcant polarization was detected in either 
22-GHz or 43-GHz observations, was used instead. The polarization position angle was calibrated by use of the 
3C286 data, applying corrections to the R— L phase difference to rotate the 3C286 polarization position angle to 
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Figure 1. Estimates of the Stokes parameters Q and U for the mapped sources, in mjy, for the two metliods described in the text 
(fitting to tile images and to tlic («, (') data). In each case, Q, in mJy, is the abscissa and U, in m,ly, is tlic ordinate, and lines on the 
plot join the {Q, U) coordinates of the first measurement to those of the second. Note that for most sources the methods agree well and 
the lines are relatively short. In cases of disagreement we provide a note in Table 2 explaining the choice of polarized flux densities used 
in the data tables. 

33°. The baseline-to- baseline scatter in this quantity was 0.6, 3.3 and 8.5 degrees for 8 GHz, 22 GHz and 43 GHz 
respectively. 

Imaging was performed using the AIPS task IMAGR with uniform weighting, on a default 128x128 grid and 
cellsizes of I'.'S, (X'? and (X'SS for 8 GHz, 22 GHz and 43 GHz respectively. For strong sources (roughly speaking, those 
containing sufficient flux density to have been used as phase calibrators) , one iteration of phase-only self-calibration 
and one iteration of amplitude self-calibration was applied. Final total intensity images were produced using 1000 
CLEAN iterations (Hogbom 1974); images in Stokes Q and U were produced using a limited number (40) of iterations 
to avoid reduction in the signal due to clean bias (e.g. Condon et al. 1998). 

Estimation of the Stokes /, Q and U flux densities was performed in two ways. The first method consisted of 
fitting to the calibrated [u, v) data, constraining the fit by fixing the position to the peak of the / map, constraining 
the intrinsic size of the peak component to be < 1", but allowing the Q and U flux density to vary. Constraints arc 
necessary to avoid instability in the fits due to noise in the polarization data. The alternative method consisted of 
finding, and fitting to, the highest absolute value in the Q and U images, within a square of side 10 pixels. Neither 
method is ideal. Fitting to the {u, v) data has the advantage of independence from effects due to deconvolution, 
notably the effects of clean bias. On the other hand, stabilising the fit requires fixing the peak of the Stokes Q and 
U to the peak of the / data, and this may not be appropriate for sources with extended polarized structure. Fig. 1 
shows a comparison of the two methods, and demonstrates that for the majority of sources the two methods agree 
well. Approximately 10% of sources show significant disagreement, mostly those with low polarized flux densities. 
In Table 2, where the data are presented, detailed comments are given in cases where the two methods disagree. 
Errors given in Table 2 are a quadrature combination of the random error, as output by AIPS, and the systematic 
calibration error, assumed to be 5% for 43 GHz and 22 GHz and 2% at 8.4 GHz. 

Because the total polarization is a quadratic sum of Q and U , it is positive-biased; values of polarized flux 
density of three times the rms noise cannot therefore be regarded as significant detections. Many studies of this 
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Frequency 1.4 8.4 22 43 

(GHz) 1.4-8.4 8.4-22 22-43 



Polarization (%) 2.2±0.2 3.1±0.2 2.7±0.2 3.1±0.2 

(1.6) (2.6) (2.2) (2.6) 

Intensity spectral index 0.01±0.04 -0.20±0.05 -0.49±0.04 

(0.04) (-0.15) (-0.37) 

Polarization spectral index 0.35±0.07 O.OliO.lO -0.09±0.06 

(0.28) (-0.15) (-0.14) 



Table 1. Average intensity and polarization quantities for the sample, counting polarization limits as zero. 1.4-GHz observations are 
from NVSS, 8.4-GHz observations from CLASS and this work, and the remainder from this work. Median values are given in brackets. 
Spectral indices are defined as = u". 

effect have been made. For example, Simmons & Stewart (1985) give a relation of estimated signal-to-noise, Po, and 
observed signal-to-noise p, in polarization for four different indicators. These all converge towards equality for p > A 
as Po = (p^ — 1)^^^, and diverge strongly for p < 2, reaching values of 1 < p < 1.5 for po = 0. The issue in these 
observations is complicated by the different values of rms noise in different images. Inspection of the images suggests 
that polarized flux densities of lOmJy are likely to be significant at all frequencies in both methods. 

Fig. 2 shows the complete total intensity and Stokes maps, together with the spectrum in total intensity and 
polarization, for the first source in order of right ascension, WMAP 0006—063. The other images can be found in 
the online material associated with this article. 

3 COMPARISON WITH OTHER OBSERVATIONS 

3.1 Consistency of polarization with the VLA caHbrator programme 

A few polarization calibrators are observed reg ularly by the VLA (Taylor & Myers 200ojl|. We compare the total 
intensity and polarization measurements of these sources with our measurements for 22 GHz and 43 GHz in Fig. 
3. Sources included in this comparison are those which were observed by the VLA in the D-configuration of 2008, 
interpolated if necessary to the date of our observation. Those sources which appear non-variable to 15% are plotted 
as separate symbols; it appears that our measurements compare well with the VLA calibration programme. 

3.2 Comparison with WMAP fluxes 

Ricci et al. (2006) comment that their AT20G flux densities from ATCA data on southern sources are systematically 
different from those of WMAP, in that the WMAP 30 GHz observations of flat spectrum sources seem to be low 
by a factor of about 1.2. Differences in this sense cannot be the result of resolution effects, since the larger beam 
of WMAP should result in a higher detected flux density in any resolved sources. We can investigate this with the 
northern-hemisphere 22-GHz data presented here, and the results are shown in Fig. 4. 

It is clear that there is a population of resolved sources, for which the WMAP flux densities at 22 and 43 GHz 
are higher than those of this work. A population of sources unresolved with either set of observations cluster around 
the line of equal flux density, although there will be some scatter induced by variability of sources between the 
two epochs of observation. At the 20% level, there is no evidence for the flux scales being significantly different 
at 22 GHz. There is marginal evidence at 43 GHz for slightly lower flux densities in these observations than in 
the WMAP observations, which may be an indication that even some relatively compact sources are beginning to 
be resolved at this frequency. It may also be a consequence of variability bias, in which some sources which are 
on average just below the 1-Jy limit, happened to be in a relatively high state when observed with WMAP. More 



Online data from the VLA polarization calibration programme is available on I http:// www.vla.nrao.edu/astro/calib/polar, 
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Image sizes: 8.4GHz: 230 arcsec, 22GHz: 89.6 arcsec, 43GHz: 44.8 arcsec. 
Flux scales; 8.4GHz (IQU) - (-2.3:2096.4, -1.0:13.3, -15.2:1.0) mjy 

22GHz (IQU) = (-4.8:1809.5, -6.3:45.8, -5.2:47.0) mJy: 43GHz (IQU) = (-9.3:1455.3, -6.9:44.5, -7.6:20.3) mJy 



Figure 2. Full data for the first WMAP source, 0006—063. The spectrum (top) is given in units of log (flux density/mjy) against log 
(frequency/GHz); error bars are shown but are smaller than the plot symbols in many cases. Filled circles are total intensity points, and 
open circles are polarized flux density points measured by fltting to the (u, v) data. Stars are polarized flux density points measured by 
maxima from the Q, U maps. The maps of the source (bottom) include, in columns from left to right, a plot of the correlated amplitude 
as a function of baseline length, maps in Stokes I, Q and U , and a plot of the polarized flux vectors superimposed upon contours of 
Stokes /. Rows, from top to bottom, correspond to maps at 8, 22 and 43 GHz. 

generally, the WMAP flux densities are significantly affected by CMB fluctuations, since a ~ 50/iK fluctuation in a 
30' beam corresponds to a flux density of about 100 mJy between 20-30 GHz. 

3.3 Comparison with CLASS 

We now compare the 8-GHz data from this project with the polarization measurements made in CLASS (Jackson et 
al. 2008), typically between 1990 and 1992, which are shown in Fig. 5 for sources with significant (6cr) polarization 
detections in both works. As expected, the total flux densities cluster around equal values for these point sources. 
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Figure 3. Comparison of intensity and polarized flux density measurements between this work and the VLA polarization calibrator 
monitoring programme (Taylor & Myers 2000) . Sources are plotted if they were observed in the D-configuration of 2008 by the monitoring 
programme, interpolated to the date of our observation if necessary. Encircled points represent sources which were stable during the 
period 3C286, 3C48, 1310+323, 1924-292, 2136+006). Dashed lines represent variation of 20% from the Une of equal flux density between 
these observations and VLA calibration monitoring observations. 

but the polarized flux densities are lower for CLASS. Again, this is likely to be a resolution effect, since the CLASS 
observations were taken in A-configuration of the VLA, which has a maximum baseline a factor of 35 greater than 
that of the D-configuration used for these observations. We pick up more polarized fiux in these D-configuration 

observations. The implication is that there is polarized emission present on arcsecond scales and therefore resolved 
out by CLASS. The relative degrees of resolution of the total intensity and polarized flux indicates that the structure 
being resolved out is more highly polarized than that which is not. 

3.4 Position angle differences and multiple components 

In Fig. 6 we plot the histograms of differences in position angle between 8-GHz and 22-GHz polarization, and between 
22-GHz and 43-GHz polarization, for objects with significant detections of polarization, have contemporaneous 
polarization measurements at all three frequencies, and have no obviously heavily resolved structure. It is obvious 
that the two lower frequencies are less well correlated in polarization position angle. Correlation is, unsurprisingly, 
slightly worse if the CLASS 8-GHz observations are used in addition to the contemporaneous 8-GHz observations. 
There is no obvious difference in position angle correlation between brighter and fainter sources. 

The obvious interpretation of these data are that significant Faraday rotation is present at the lower frequencies; 
since this effect is proportional to A^, the angle of rotation would be expected to be about 6-7 times larger between the 
lower two frequencies. Galactic Faraday rotation is unlikely to be responsible for this effect, since the <100 radm~^ 
typical of the Galactic foreground is an order of magnitude smaller than the rotations needed significantly to affect 
8-GHz measurement, and indeed in our data there is no correlation of this rotation with galactic latitude. Zavala & 
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Figure 4. Comparison of the 22-GHz and 43-GHz flux densities in this work with those of Wright et ah (2009; quoted to two significant 
figures in that paper). Lines are drawn at 20% either side of the line of equal flux. 
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Figure 5. Comparison between flux densities in this worlc and CLASS of total fiux density (filled circles) and polarized flux density 
(open circles) for cases where both surveys detect significant polarized flux. Note the generally lower polarized flux density in CLASS. 

Taylor (2004), in a study of parsec-scale polarization of quasar radio cores, find core RMs from 500 to a few thousand 
rad m~^, and suggest that these are due to a foreground Faraday screen close to the radio jet. Such parsec-scale 
components are likely to dominate the 22- and 43-GHz flux densities. On the larger scales probed by the lower- 
frequency observations, it is possible that in many sources there are multiple polarized components with different 
position angles. The fact that more polarization is detected here in D-configuration than in the A-configuration 
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Figure 6. Differences in position angle between observed polarization at different frequencies. Solid line: difference between 22 GHz 
and 43 GHz. Dashed line: difference between 8 and 22 GHz for contemporaneous observations. Dash-dotted line: difference between 8 
and 22 GHz for objects where only archival CLASS 8-GHz data is available. Only observations with significant detections of polarized 
flux (defined as 6(T in CLASS and as objects with visible polarized flux in the Q and U maps in these observations) are plotted. Known 
heavily resolved sources are excluded. 



CLASS observations suggests that a more extended jet component, with a lower RM than the core, may be present. 
Future high-sensitivity observations with e-MERLIN and the EVLA should reveal any such low surface-brightness 
components. 



4 CONCLUSIONS 

We have presented observations at 22 and 43 GHz of 199 of the 203 radio sources in the WMAP5 list of Wright et al. 
(2009) which have flux density >lJy at 22 GHz, and observations at 8.4 GHz of 133 of these objects. The total flux 
densities are broadly consistent with those listed by WMAP. For the few objects with previous VLA observations of 
polarized flux densities, there is also good agreement. There is evidence for extra polarized flux density on arcsecond 
scales at 8.4 GHz compared to earlier observations at higher resolution. Polarization position angles are very similar 
at 22 and 43 GHz, but are much less well correlated with position angles measured at 8.4 GHz. This is likely to 
be due to multiple polarized components which make up the generally complex radio spectra of flat-spectrum radio 
sources. In a further paper we discuss correlations between the measured intensity and polarization parameters in 
more detail, and calculate the implications for future CMB polarization measurements of the distribution of polarized 
flux densities in these discrete sources. 



ONLINE MATERIAL 



The complete maps can be accessed at http : //www . j b . mein . ac . uk/ ~nj j /pol jnaps/l This figure includes, for each 



source, in columns from left to right, a plot of the correlated amplitude as a function of baseline length for each 
source, maps in Stokes /, Q and U, and a plot of the polarized flux vectors superimposed upon contours of Stokes 
/. For each object, rows are given corresponding to maps at 8, 22 and 43 GHz respectively. 

Plots of spectra can be accessed at http : / / www ■ j b . man . ac . uk/ ~nj j /pol_spec/| and are given in units of 
log(flux density/mJy) against log(frequency/GHz). Filled circles are total intensity points, and open circles are 
polarized flux points measured by fitting to the (u, v) data. Stars are polarized flux points measured by maxima 
from the Q, U maps. 
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Caption to Table 2. 

Flux densities of observed sources at X-, K- and Q-bands. Fluxes in italics are from CLASS, in the case of multiple observations 
using that with the lowest error in polarized flux density. A dagger indicates that the measurement has been taken from the Q and U 
maps, rather than the fits to the u-v data (see text for details). A "zero" error means smaller than 1 mjy. Notes to the table: 1. Nearest 
plausible ID is 40+25.01 (00 19 39.21 +26 02 45.4) which is 0.5 Jy at 1.4 GHz (NVSS); multiple sources may contribute to the WMAP 
flux density. The NVSS position was observed, leading to the detection of a resolved, polarized source at X-band. 2. Source is identified 
with 3C20 which is a 50" steep-spectrum double (Laing 1981). Nearly all structure is resolved out by these observations. 3. Identification 
is the star-forming galaxy NGC253. Structure is nearly resolved out by these observations. 4. No X-band polarization available from 
CLASS. 5. Identification is 3C33 which is a 5' steep spectrum double (e.g. Leahy & Pcrlcy 1991). Only the central regions arc seen in 
these observations. 6. Used as the zero-polarization calibrator at X-band. 7. Offset polarized component? Fit to the K-band images gives 
13.4 mJy polarized flux density at 11°. No significant difference between methods at Q-band. 8. CLASS X-band polarization measured 
on two occasions: 21 mJy in PA 57 (total flux density 1039 mJy), 29mjy in PA 25 (total flux density 1061 m.ly). 9. Identification is the 
pair of sources 3C 66A, 3C 66B. Observation pointed at 3C 66A and structure heavily resolved; most of the high-frequency contribution 
to the flux is probably from 3C 66B, however. 10. No X-band polarization available. K-band polarization from fltting to ti — ?; data 
probably underestimated; images give 10.2 mJy in PA 25°. 11. Source is identifled with 3C78 (e.g. Saikia et al. 1986); heavily resolved 
by these observations. 12. Identifled with 3C84; used as zero-polarization calibrator for K-band and Q-band observations. Polarization 
level difficult to measure because of residual errors; best guess from analysis of images is 13 mJy polarized flux density in PA 129° 

13. Identified with 3C98 which is a 216 arcscc steep spectrum double (e.g. Leahy et al. 1997); heavily resolved by these observations. 

14. Probably resolved; peak in the Q-band image has a polarized flux density of llmjy at PA 0°. 15. Extreme discrepancy between 
flux observed in all VLA bands and the WMAP flux. No obvious identiflcation from NVSS; the observed position (04 24 2.58 +02 26 
42) is the closest NVSS source (about 10') from the WMAP position; WMAP flux may be a combination of this and two other ~1-Jy 
sources at only slightly larger distance. 16. Fits to K-band image suggest 12mjy polarized flux density in PA 117°, but sparse data 
make estimation difficult. 17. Blended source (see note in Wright et al. 2008). 18. ID is a planetary nebula (IC0418). 19. Several X-band 
measurements from CLASS with highly variable polarized fiux density and PA. 20. Identification is a very large (15') double-lobed radio 
source, PKS0634— 20. One of the lobes was observed, and is heavily resolved. 21. Possible offset polarized component in K-band; analysis 
of image suggests lOmJy in PA 14° . 22. ID is probably the strong source 4.5' SSE, 3C196. This is a small double (e.g. Reid et al 1995) 
both of whose components are detected. 23. Very little data obtained at X-band and K-band, no reliable results at these frequencies. 24. 
Phase errors limit reliability; analysis of images suggests K-band polarization 44m Jy at PA 23°. No significant polarization detected at 
Q-band. 25. Extremely steep spectral index between K-band and Q-band. 26. ID is Hydra A, a large radio source (e.g. Taylor et al. 1990). 
Core of the source was observed, consisting of extended feature 45" across. Polarization level depends strongly on aperture used. 27. 
Analysis of images suggests higher K-band polarization (50 mJy) at same PA. 28. Identification is 4C39.25. X-band polarization (CLASS) 
highly variable. 29. Identification is M82 (starburst galaxy). 30. Error in pointing of array. 31. Clear ID, but data of poor quality, phase 
solutions not good. Polarization measurements unreliable at K and Q-band. 32. Poor quality data at K-band and Q-band, phase solutions 
not good. Polarization measurements unreliable at both frequencies. 33. Most plausible ID is 4' W, little data at K-band and phase 
solutions not good. 34. Good ID, poor quality data at K-band and Q-band, phase solutions not good and polarization measurements 
unreliable. 35. Good ID, poor quality data at Q-band, phase solutions not good and polarization measurements unreliable. 36. Analysis 
of Q-band image suggest slightly lower polarized flux density, 12m Jy at PA 80°. 37. ID is 3C270; pointing position is that of the core 
(Birkinshaw & Davies 1985). 38. ID is 3C273; polarization at X-band (CLASS) is high and variable. 39. ID is M87; the polarized jet 
pointing WNW is visible. Very little data at K-band; this polarization measurement is unreliable. 40. Sparse data at K and Q-band; 
image at Q-band suggests that this polarization measurement is unreliable. 41. Fits to Q-band images suggest polarized fiux density of 
12mjy in PA 61°. 42. Little data at K-band; unreliable polarization measurement at this frequency. 43. ID is 3C286; primary flux and 
polarization position angle calibrator; also used for non-closing offset calibration at X-band. 44. IC4296. This is 2.5 arcmin double with 
compact core. No significant polarization detected. 45. Possible offset polarized component in K-band; analysis of image suggests 14mjy 
at PA 66°. 46. No obvious catalogued source. Nearest radio source is MRC 1424-328 (NVSS position 14 27 35.7, -33 02 22, 253 mJy 
at 1.4 GHz) approximately 1.3 arcminutes E. This position yields 180 mJy of correlated flux at 8.4 GHz and no detections at higher 
frequencies. WMAP flux density of 1 Jy at 22 GHz may be the sum of a number of nearby sources. 47. Analysis of K-band image suggests 
slightly higher polarization: 13mjy at PA 75°; detection is significant. 48. ID is NGC6251; position of core observed (e.g. Perley, Bridle 
&c Willis 1984). 49. ID is 3C345. Phase errors in K-band observation; analysis of image suggests slightly higher polarized flux density 
of 75mjy in PA 164°. 50. ID is Hercules A; core observed (e.g. Gizani & Leahy 2003), highly resolved. The core lies in between two 
polarized components, and its polarized flux density is much less than that inferred from the u — v data by the automatic fltting. 51. 
Fits to Q-band images suggest polarized flux density 14mJy in PA 60° . 52. ID is 3C380. 53. K-band images reveal significant detection 
of polarized flux density, ~10mjy in PA 117°. 54. Possibly offset polarized component at K-band of 16mjy in PA 62°. 55. ID is 3C454.3. 
56. K-band polarization measurement unreliable; noisy map with signiflcant phase errors. 



